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The question of the origin of fermion mass hierarchy, mixing angles and CP violating 
phase is one of the most outstanding problems in particle physics. There have been numerous 
attempts to study this problem, some of which are more theoretical in nature while others 
are more phenomenological. However, it is generally agreed that the final word is far from 
being said. Furthermore, it is also agreed that the solution, whatever it might be, is to be 
found outside of the Standard Model (SM). 

In all of these studies, the phenomenological-ansatz approach is much more modest in 
scope. Starting with some simple assumption about the form of the mass matrix whose 
theoretical justification is yet-to-be-determined, one could fit quark masses and mixing 
angles. One of such approaches is particularly appealing: The pure phase mass matrix 
(PPMM) 0] . This particular ansatz is based on a simple assumption that there is a single 
and unique Yukawa coupling for each quark sector and that the 3x3 mass matrix takes the 
form Ai = gy^v/ ^/2){exp(i6ij)}, where i,j = 1,2,3. This kind of mass matrices belongs to 
a class of the so-called "democratic mass matrices" (DMM) [Q. The pure phase mass matrix 
is attractive in that the hierarchy of masses is governed by a single Yukawa coupling in the 
limit where all phases vanish. A realistic hierarchy comes about when the phases, which are 
treated as small perturbations, are put back in. Although it is conceptually attractive, no 
attempt was made to justify its underlying assumption. Earlier works on trying to model 
the pure phase mass matrix relied entirely on the framework of four-dimensional field the- 
ories. Although there are a number of useful lessons that can be learned from this mode of 
thinking, one is sometimes faced with more questions than answers. 

On another front, there has been important conceptual developments in the last few years 
related to a possible existence of Large Extra Dimensions §]. Not only does this concept 
force us to rethink about notions such as the question of what the ultimate fundamental scale 
of nature might be, it also inspires us to reformulate some of the longstanding problems in 
particle physics such as the origin of fermion masses and mixings. The hierarchy of masses 
has been reexamined recently within the framework of large extra dimensions, and new 
interesting ideas have emerged such as the notion of "thick branes" and the localization of 
various fermions inside these branes . This localization can be accomplished by a domain 
wall inside the brane. This gave rise to the idea of the strength of the Yukawa coupling 
(which is proportional to the mass of the fermion) as being the overlap of the wave functions 
of the localized fermions. As stated in Ref. it is easy to think of the reason why some 
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fermions are heavy and some are light: The heavy ones have large overlap and the light ones 
have small overlaps. There has been some works done along that line in order to explain 
the fermion mass hierarchies. Most of these works made use of the size of the wave function 
overlaps to discuss the fermion mass problem. 

Whatever various scenarios might be, the common important elements which transpired 
from these works are basically the locations of the domain walls and the size of the wave 
function overlaps. In fact, many of the physics results will depend on the actual placements 
of the domain walls along the extra dimensions. 

Our approach in this paper is as follows: For each fermion sector (e.g. the up and down 
quark sectors), there is a universal overall mass scale whose Yukawa coupling strength is 
determined by the size of the overlap. This gives rise to a democratic mass matrix whose 
elements are all equal to unity, apart from a common mass scale factor multiplied by an 
effective Yukawa couphng. All that is needed is to locahze all the left-handed fermions at one 
location, regardless of family indices, and all the right-handed fermions at another location 
along the fifth dimension inside the thick brane, and, in addition, to endow the fermions 
with a permutation symmetry. Unfortunately, it is well known that this kind of matrix does 
not work: one obtains one non-zero mass eigenvalue and two zero eigenvalues. The matrix 
{1} has to be replaced by another quasi-democratic one of the form such as {exp{i9ij)} 
for example. The mass hierarchy which arises within each sector is due, in our scenario, 
to the introduction of a sixth dimension and a thick brane along it. The introduction of 
"family" domain walls at different locations inside this thick brane generate different phases 
for different famihes. It will be seen that it is these phase differences which give rise to the 
pure phase mass matrix or, as we shall see, an almost-pure phase mass matrix. 

We would like to make the following remark. Our model will contain a certain number 
of parameters that need to be fixed phenomenologically. However, what we present here 
is a new perspective on an old problem which, hopefully, can give further insights which 
might be useful for future investigations. What we are doing here is to try to rephrase the 
origin of quark mass hierarchy (and eventually that of the leptons as well) and CP phase 
in a completely new context: that of the Compact Extra Dimensions (CED). We will show 
below that the appearance of the phases in the mass matrices, a crucial element in their 
construction, appear rather "naturally" . From this point of view, it appears to be a definite 
conceptual advantage of the CED scenario. 
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One remark is in order here concerning the introduction of a sixth dimension. It is well 
known that, with just one extra compact dimension, the fundamental 5-dimensional Planck 
scale cannot be of the order of a few TeV or so, for it will introduce deviations to the inverse 
square law on astronomical distances. Recent gravity experiments ^ down to a millimeter 
or so put a lower bound of around 3 TeV on the 4+n Planck scale for the case of n=2 
(with equal compactification radii). This fact, of course, was not the one motivating us 
in introducing a sixth dimension. It is rather the natural way in which phase differences 
appear between different fermions eventually giving rise to a pure phase mass matrix which 
motivated us. 

The organization of the paper is as follows. First we review various features of fermions 
in five dimensions, including, for instance, the concept of fermion localization. We then 
show how, with a rather simple assumption, a democratic mass matrix appears. Next, we 
introduce fermions in six dimensions and show how phase differences appear, and how one 
can construct an (almost) pure phase mass matrix from this result. In this construction, 
"family" domain walls are introduced and it is shown that their small separations along 
the sixth dimension are responsible for the aforementioned phase differences. Unlike what 
happens along the fifth dimension, the fermion wave functions are not of the localizing type 
but are rather oscillating. We will then discuss how hermitian and non-hermitian pure phase 
mass matrices arise. Finally, we will discuss some possible connections to the strong CP 
problem [^. 



I. FERMIONS IN 5 DIMENSIONS AND DEMOCRATIC MASS MATRIX 



A. A Review 



In this section, we will review some aspects of fermions in five dimensions which have 
support [0, L] along the fifth dimension. In other words, we are discussing a "thick brane" 
of thickness L. This discussion serves two purposes: to set the notations and to lead to the 
democratic mass matrix. 

We will adopt the effective field theory approach of Refs. [ill . This approach has the 
merit of being relatively simple and transparent as far as the physics is concerned. We first 
summarize below what has been done for the case of one flavor of fermions, without and 
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with a background scalar field. 

To set the notations straight, the 4- dimensional coordinates will be labeled by x^^ with 
/i = 0, ..,3 while the fifth coordinate will be labeled by y. We start out with a free Dirac 
spinor of 5*0(4, 1) which has four components, tp. The gamma matrices are and 7^ = ^75. 
The free Dirac Lagrangian is given by 

The above Lagrangian has the following Z2 symmetry: tp{x, y) —>■ "^{x, y) = ±75-?/'(a;, L — 
y). When this symmetry is combined with the periodic boundary condition: ip{x^y) = 
\E'(x, L + y) = ipi^x, 2L + y), one obtains: ip{x, —y) = \E'(a;, L — y) = ±j^iIj{x, y) and ip{x, L + 
y) = '^{x,y) = ±75'?/'(a;, L — y), which shows that y = 0,L are fixed points. One can 
subsequently define the chiral components of ip by using the usual operators Pr^l = (1 ± 
75)/2, with ip^ = Pj^ip and ip^ = Pl'4'i with 75'?/'^ = iV'it- The previous symmetry and 
boundary conditions are what usually referred to in the literature as compactification on 
an S1/Z2 orbifold. One can have fermions which have the symmetry il){x,y) —>■ \E'(a;,2/) = 
+75'?/^(x, L — y), and those which have ipi^, v) ~^ "^{x, y) = —■y^ipi^x, L — y). 

For simphcity, we shall discuss the case ip{x, y) y) = +75'?/'(a;, L — y) below. This 

corresponds to the case where only right-handed zero modes survive in the brane, as shown 
below. For the other situation, ip{x,y) — > \l'(x,?/) = — 75'?/'(x,L — y), only the left-handed 
zero modes survive inside the brane, as one can easily check. 

Zero modes residing in the brane are supposed to be independent of the extra coordinate, 
y in this case. From the above discussion, one can see that ip- vanishes at the fixed points, 
and hence there is no zero mode for ip^. The only non- vanishing zero mode is ipo+- This 
can also be seen explicitly by writing 

ipM+ ix,y) = ipM+ ix)^M+ (y) , (2a) 



ipM-{x,y) = ipM-{x)^M-{y), 



(2b) 



for a mode of mass M. From the explicit solutions for ^ as given in Ref. [0], one can again 



see that there is only one chiral zero mode inside the brane. Four-dimensional chirality is 
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seen to arise from the symmetry and boundary conditions. The chiral zero mode il'o+ is 
uniformly spread over the fifth dimension y. To locahze ?/'o+ at specific points along y inside 
the brane, the use of domain walls have been suggested by Refs. 0, To this end, a 
background scalar field, $, is introduced. The Lagrangian is given by 

£ = ^ (^^ ^ - 75^ - + (3) 

The symmetry and boundary conditions on $ are now: $ — > L — y) = — y); 

—y) = L — y) = —^{x, y) and L + y) = $(2;, y) = — L — y). It can then 
be seen that vanishes at the orbifold fixed points: y = 0,L. As discussed in Ref. $ 
has a minimum energy configuration: y)) = (f){y), with 0(0) = (p^L) = 0. From the 

modified equations for ^m± with an added term f4>{y), one can easily see the localization of 
the zero mode, namely 

^,4y)=ke-'(y\^,_{y)=0, (4) 

where 

s{y) = f I'dyW). (5) 
^0 



As pointed out by Ref. [T^, the chiral zero mode, ^o+{y), is now localized either at y = Q 
or y = L depending on the sign of f(f){y)- 

As in Ref. 0, the special choice f(f){y) = 2fi'^y which makes the operators a = dy + f(f>{y) 
and a"!" = —dy + f4'{y) behave like the annihilation and creation operators of a Simple 
Harmonic Oscillator (SHO), the normalized wave function for the chiral zero mode C,o+{y) 
takes on the famihar form C,o+{y) = (y^/(7r/2)^/'') exp(— /i^y^). One clearly notices the 
localization of C,o+{y) at y = 0. Another way of describing this phenomenon is the fact that 
has a kink solution of the form V taioh({\/2)^^'^Vy) which basically traps the fermion to 
a domain wall of size {{X/2y/^V)-^ p[. 

The next question concerns the possibility of localizing the chiral zero mode at some 
other location than the one at the orbifold fixed points. Ref. has proposed to change 
the Yukawa interaction '0(/0(t/))-0 to ip{f(f){y) —rnjip so that the wave function of the chiral 
fermion field is now localized at the zero of f(f){y) — rn instead of /0(?/). With the SHO 
approximation, this zero would he at y = m/2/i^. However, in order to be compatible with 
the Z2 symmetry of the Lagrangian, as shown in Eq. (R[), one should also require a "mass 



reversal" m —m simultaneously with the Z2 transformations. This is the assumption we 



will be making in this manuscript. (Another approach is given in Ref. |T0|). 

As emphasized by Ref. different massless chiral fermions can be localized on dif- 
ferent slices along y, inside the thick brane. These locations are determined by the ze- 
ros of f(j) — rrii = 0. Within the SHO approximation, the wave functions are given by 
{y/JI/ {n/2Y^^) exp(— — yi^), where yi = mj/2/i^. The interesting idea proposed in Ref. 
IP is that the effective Yukawa couplings between SM fermions and SM Higgs scalar, which 
eventually determines the size of the mass term, are mainly determined by the wave func- 
tion overlap between the left-handed and right-handed fermions. Hierarchy of masses then 
appears to depend on the size of the overlaps. 

From hereon, we shall turn our attention to left-handed zero modes inside the brane as 
used in the SM. As we have mentioned earlier, these come from five-dimensional fermions 
with the Z2 symmetry ipi^, y) y) = — 75'?/'(x, L — y). 

To prepare the groundwork for our subsequent discussion, let us write down the action 
in five dimensions of a left-handed fermion, a right-handed fermion, and the Yukawa in- 
teractions with a background scalar field, and a SM Higgs field. Following Ref. |0, we 
will denote quarks in five dimensions by the five-dimensional Dirac fields: {Q, If^, D'^) and 
their left-handed zero modes by the following Weyl fields: {q,u'^,d'^). Notice that with this 
notation, a right-handed down quark, for example, will be d'^. Since we will be dealing in 
this paper solely with the quark sector, we are not writing down the lepton fields. This 
will be dealt with in a subsequent paper. The SM transformations of the above fields are 
self-evident by the use of these notations. In addition, one also introduces two sets of scalar 
fields: a SM singlet background scalar field, 0, whose VEV is y)) = 4>{y), a SM doublet 
Higgs field H{x, y) whose zero mode h{x) is assumed to be uniformly spread along y inside 
the thick brane. The 5-dimensional action can be written as 



S = j d'xQit h + ]m)Q + h + - ^vW (6) 

where C5 = joj2jy From the above equation, one notices that Q, U^, are localized at 
Uq = ^1 yu = '"^c//2/i^, yo = mD/^fi"^ respectively. In principle, mu and ku can be different 
from and kd respectively. However, as we can see below, it is sufficient to have mu ^ mjj 
in order for the resulting masses of up and down quarks to be different, even if ku = kd- 



Assuming that the zero mode of H is uniformly spread over y inside the thick brane, the 
4-dimensional effective action for the Yukawa interaction for the up quark can be written as 



and similarly for the down quark. From the form of the wave functions, one obtains the 
4-dimensional effective Yukawa couplings for up and down quarks as follows 



Two remarks can be made concerning Eqs. (|]) and (j^). First of all, as emphasized by 
Ref. 0, even if k's are of order unity, the effective Yukawa couplings can be quite small if 
fJ'yu,D ^ 1- Basically, the size of the effective coupling is sensitive to the relative distance 
between left and right-handed quarks as compared with the characteristic thickness of the 
domain walls. The second remark concerns the Yukawa couplings in five dimensions, ku^d. 
In this new framework of large extra dimensions, one has to separate the mechanism which 
separates qy^u from gy^d, already at the level of the 5-dimensional action from that which 
separates gY,u from gy^d at an effective field theory level in four dimensions due to different 
localization points along the extra dimension inside the thick brane. It might happen that 
the 5-dimensional action has an up-down symmetry in the Yukawa sector which is broken 
down inside the brane. We shall return to this question at the end of the paper. 

B. Democratic Mass Matrix 

Let us, for now, concentrate on just one sector, e.g. the up sector. Let us assume that 
there are three families. The fermion fields in five dimensions that we will be dealing with in 
this section will be Q and f/^. As we shall see below, in order to obtain the DMM scenario, 
we will put all the Q's at one location along y inside the thick brane, and all the f/'^'s at 
another location. With this simple assumption and the assumption that the SM Higgs zero 
mode is uniformly spread inside the thick brane, one can naively obtain the democratic 
mass matrix mentioned above. However, with the gauge field zero modes also spreading 






(8) 




(9) 
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uniformly inside the thick brane, this will give rise to unwanted flavor- changing neutral 
current (FCNC) operators. A symmetry has to be imposed in order to avoid these FCNCs. 

A simple symmetry that one can use is a permutation symmetry among the three families, 
for both Q and f/^ One can have: O ^3^', with Q S^Q and f/^ S^'U^. The 
background scalar fleld described earlier (j){y) is a singlet under the above permutation 
group. (In this way, one will see that all Q's are localized at one place and all f/'^'s are 
localized at another place.) One can now include gauge interactions in the kinetic terms of 
(^) by making the replacement ^5 — ^Ps, namely 

So = j (fx Q{i P, + f<l){y))Q + U-{t A + f<P{y) - mu)U' (10) 

It is simple to see that Sq is invariant under the above permutation symmetry. Eq. (p!OD also 
implies that all Q's are localized at one place and all If^'s are localized at another place. 

Next, we wish to introduce a Yukawa interaction between the SM Higgs scalar and Q 
and If^. First, we notice that a term such as 

C-yukawa = HuQ^C^HW^ + h.C. (11) 

breaks the permutation symmetry since Q and transform under different groups. If they 



were to transform under the same permutation group, Eq. (]TT|) would be an invariant. 
However, it would give a mass matrix of the form 

'^l 0^ 



1 



(12) 



which is not of the DMM type. It turns out that with S'^ ® S]^" , one can construct an 
invariant for each permutation group: '^iQi for 5*^, and for 5*^" where i = 1,2,3 

and j = 1,2,3 are family indices. From this, one can construct an invariant action for the 
Yukawa interaction 



Sii 



^yukawa 



The effective action in four dimensions can now be written as 



j (fx Ku QjC^H J2 Uj + h.C. (13) 



Seff, Yukawa = J <fx Ku ^ q'^'\x)h{x)u'''^ j dyCq{y)ii<^{y) + h.C. (14) 



Since all the g^'s are located at the same place inside the brane, and similarly for all the wf, 
the wave function overlap / dy ^^g{y)^l^c{y) is universal and independent of With this, 
one can now rewrite Eq. (|T3) as 



Seff, Yukawa — \ d'^X Qy^uQ^ [x] 



1 1 1 
1 1 1 



h{x)u'' + h.c. 



(15) 



where gy^u is given by Eq. (|), = {qf, gj, q'^) and similarly for u^{x). From Eq. one 
obtains the democratic mass matrix 

/i 1 1\ 

1 1 1 



111 



(16) 



An important remark is in order here. The universal strength in Eq. (0) depends on, 
besides the SM quantity v/\^ ~ 175 GeV, gy^u which is a product of two factors: the five- 
dimensional Yukawa coupling, k, and the overlap of left-handed and right-handed fermion 
wave functions. In this scenario and its extension presented below, it is this product that is 
important, and not simply the size of the overlap. 

As we have mentioned above, the above matrix can be brought by a similarity transfor- 
mation to a form 



M' = SMS-^ 



(17) 



'^O 0^ 



9y,u 



V2 
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As one can see above, one needs to move beyond the DMM scenario in order to obtain a 
more "realistic" mass matrix. This is what we propose to do in the next section. 

One might wonder what the distinctive feature a fifth dimension has to give us in regards 
with the above problem. Could one not obtain a similar result staying in just four dimen- 
sions? In principle, the answer is yes. However, it appears more attractive to think that, 
once qi are lumped together at one place and are lumped at another place, one would 
obtain the DMM naturally. It is interesting to envision a scenario in which the Yukawa 
couplings are as universal as the gauge couplings themselves, with the possibility that the 
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effective Yukawa couplings can be different from one another due to the different overlaps 
between left and right fermions. (Gauge interactions are chirality conserving and, as a re- 
sult, the effective gauge coupling with the gauge boson zero mode is the same as the original 
couphng.) 

The above discussion carries over to the down sector in a similar fashion. Obviously, 
although attractive, this kind of democratic mass matrix does not give the correct mass 
spectrum. An extension of DMM was discussed by Ref. in which, instead of having 
one's as matrix elements, one has pure phase factors such as exp{i6ij). (The diagonal 
elements can be all unity by a suitable redefinition of the quark phases.) Explicitly, a pure 
phase mass matrix looks like = gy^v / \/2){exp{i6ij)). 

To construct a model for PPMM- even for the special case such as a symmetric matrix, 
one usually requires a rather complicated Higgs structure 0. That is if one stays in four 
dimensions. One might wonder if extra dimensions might help in this regards. We have seen 
above how an additional dimension could help conceptually in obtaining a democratic mass 
matrix. The question we ask is the following: Could pure phases such as exp{i6ij) arise from 
extra dimensions and not from some kind of complicated Higgs sector? In particular, if we 
keep the Higgs sector to a minimum (one Higgs), this phase cannot come from the Yukawa 
coupling nor from the VEV of the SM Higgs. We have seen that, in five dimensions, a chiral 
zero mode has, as a part of its wave function, ^{y) which behaves, upon being trapped by a 
domain wall, like exp(— /i^y^). As we shall see below, by adding another compact dimension 
(the sixth one), the phases appear as the overlaps between wave functions of fermions which 
are "trapped" at different locations along the 6th dimension. What this really means will 
be explored in the next section. 

II. FERMIONS IN 6 DIMENSIONS AND PURE PHASE MASS MATRIX 

Notwithstanding the string theory argument, there might be another simpler motivation 
for the need of more than one extra spatial dimension: If the fundamental 4 + n "Planck" 
scale were of O(TeV) to "solve" the hierarchy problem, and if the n extra dimensions were to 
be compactified with the same radius R then n > 2 in order for R to be in the submillimeter 
region as required by the lack of deviation from the ordinary inverse square law down to 
about 0.2 mm [0. In our case, the above need is dictated by the desire to build a more 
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"realistic" mass matrix: the so-called pure phase mass matrix or its almost-pure-phase 
counterpart. (In this construction, we are not concerned about whether or not the ultimate 
theory contains more than six dimensions.) To this end, we first study the behaviour of 
fermions in six dimensions, subject to similar boundary conditions as in the 5-dimensional 
case. 

A. Fermions in six dimensions 

The task of this section is to study fermions in six dimensions, with the ultimate aim of 
obtaining massless chiral fermions in four dimensions. 

In order to discuss fermions in six dimensions, we first turn our attention to the represen- 
tation of gamma matrices for these fermions. Before we begin the discussion, a few remarks 
concerning spinors in SO{N) are necessary. 

We shall be working with the group S0{5, 1) that, as we discuss in Appendix ^ has two 
irreducible spinor representations of dimension 4. We shall put tp+ and tp- into a reducible 
"Dirac" spinor ip = The chiral representation of the gamma matrices for S0{5, 1) 

is shown in Appendix ^ The notation for the coordinates will be similar to the five- 
dimensional case, with the sixth dimension denoted by z, namely xn = {xo,xi,X2,X3,y, z). 
The free Lagrangian for ip is now written as 

C^ = tijT^dN^, (18) 

where = 0, l,2,3,y,z. The metric used in this paper is simply (-+++++). It is useful 
to see explicitly the Lagrangian written in terms of the components of ip. For this purpose, 
we give the explicit forms for Ty and as can be seen from Appendix 0, 

r.-(' (19) 




(20) 



E 

1 

where 75 is the usual matrix encountered in four dimensions and is a 4 x 4 unit matrix. 
In addition, we also need ip = V'^To = ("^-i ~'^+)- Eq. (^) can now be rewritten as 

= -i ^lj+Ydt,i>+ - i ^/j-^y^di^i/j. + ilj+^^dy^lj+ + ilj_-f^dyi>_ - i ^lj+d^i>+ + i i>-d^^lJ- ■ (21) 
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As we explain in the Appendix, the 4-dimensional kinetic terms (the first two terms of the 
above equation) will acquire a plus sign when 7^ are replaced by 7^^ which are appropriate 
for the metric ( — h ++) which is a remnant of the original metric ( — h + + ++)• The reader 
is strongly recommended to consult Appendix ^ concurrently with this section in order to 
avoid confusion. 

As in the case of the fifth dimension, we will assume that the sixth dimension is compact- 
ified on an orbifold S1/Z2. ip is assumed to have support [0, Lg] along the sixth dimension. 
We first discuss this Z2 symmetry for free fermions. 

From Eq. (|18|), one can see that the Lagrangian has the following Z2 symmetry: 

z) z) = r^V(a;", Lq-z). (22) 

With Fj, given above, this symmetry translates into 

iJ+{x'',z) ^ ^+(x",z) = V^_(x°,L6-2;), 
iIj_{x",z) ^_(x",^) = i/j+{x'',Lq- z) . 

As with the five-dimensional case, our boundary condition is 

7/>±(a;°, z) = ^±(a;", + z) = V'±(x", 2L<, + z) . (24) 
Again, combining (|23|) with (|5), one obtains 

^±(x",-z) =V^^(a;",z), (25) 



(23) 



^±(x^L6-^) =^^(x^L6 + ^). (26) 

We immediately recognizes 2; = 0, Lg to be the fixed points of the orbifold. It is convenient 
to rewrite ip± as 

^± = ^ix±v)- (27) 
In terms of x and 77, the boundary conditions become 

X{x",-z)=x{x",z), (28) 



?7(x°, —z) = —r]{x"', z) 
13 



(29) 



(30) 



77(x°,L6-z) = -r/(x°,L6 + z). (31) 
From the above boundary conditions, one can see that r] vanishes at the fixed points z = 

As usual, we shall write: 

XMix",z) = XMix"')^^,Miz), (32a) 

?7M(a;", z) = ?7Af(a;")^^,M(-2)- (32b) 
Since the zero modes in the "4-brane" are independent of z, we have 

x(a:",^)o = A;x(x°);r7(x°,z)o = 0, (33) 

where A; is a constant. Again, the free fermion wave function for the zero mode is uni- 
formly spread over the 6-th dimension. We now investigate the effect of a coupling with a 
background scalar field having a kink solution. 

For the discussion which follows, it is convenient to notice that 

-i ip+d^ip+ + i dzip_ = i xdzV -iv^zX- (34) 

Eventually, we would like to find an equation for the surviving zero mode x{x", z)o in the 
presence of a background scalar field which will be assumed to be real. For this purpose, let 
us write the surviving zero mode x 

Xo(x",z)=x(x°)ex,o(^). (35) 



As we shall see, upon using Eq. (34) and subsequent interaction terms, one can derive an 



equation governing the behaviour of ^x.o(^) along z which will eventually tell us whether or 
not one has a localized behaviour as in the five- dimensional case or an oscillatory one (pure 
phase). This will depend on the type of fermion bilinears which couple to the background 
scalar. Roughly speaking, if the coupling ends up to be of the form if]xh{z), for example, 
then C,x,o{z) will have an exponentially-suppressed form similar to the five-dimensional case. 
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If, however, it ends up looking like f]xh{z), then Cxfii^) ^^^^ have an oscillatory behaviour. 
This is so because of the way Eq. (0) looks. 

We now look for the aforementioned fermion bilinears which are required to be hermitian 
(because the background scalar field is assumed to be real) and Lorentz invariant. 

Let us introduce a real scalar field which transforms under Z2 as 

$(a;",z) -^-<l>(x",L6-2). (36) 

First, the most obvious, hermitian and Lorentz-invariant bilinear is simply (remembering 
that Tq is anti-hermitian with our metric) 

i^{x'',z)ij{x'^,z). (37) 

Notice that (|37|) , when expanded in terms of x and rj, are of the form ifix + ■■■■ This, when 
combined with Eq. (|34D, would give an exponentially-suppressed form for the zero mode if 
there exists such a Yukawa coupling. Can it couple to $? If the reflection Z2 symmetry 
were the only symmetry around, it is straightforward to see that a coupling of the form 
i il){x'^ , z)iIj{x°' , z)^{x°' , z) is an invariant. This, as we have mentioned above, would not be 
what we are looking for, namely an oscillatory wave function. A mere mimicking of the 
five-dimensional case would not work. Below we propose a mechanism where the desired 
behaviour could arise. 

Let us endow the scalar and fermion fields with an additional discrete symmetry which 
will be called the Q-symmetry and which works as follows. Let us divide the space inside 
the brane of thickness Lq into two regions: to ^5/2 (Region I) and Lq/2 to Lq (Region II). 
Let us define the following transformations. Under Q, 

<l>(x",^) ^ -$(x",^). (38) 

Notice that ( pSf ) is not to be confused with (|36D which is a reflection symmetry. We then 
notice the following fact: If z is inside Region I then Lq — z will be inside Region II and 
vice versa. For the fermion, we will impose the following Q-transformations: ip ^ ip for z 
in Region I and ip —ip for z in Region II. 

With the above Q-symmetry, one notices that a coupling of the form 
iIj{x"' , z)ijj{x°', z)^{x°' , z) is forbidden for any point z inside the brane. However, a 
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nonlocal interaction of the form 'iplx"', z)iIj{x°' , Lq — z)^{x°',z) is allowed by the Q- 
symmetry. In particular, a hermitian bilinear containing iplx"', z)iIj{x°' , Lq — z) of the form 
z)ip{x", Lq — z) — Lq — z)ip{x°', z) is allowed by this symmetry. 
The way the Q-symmetry works seems to imply that the orbifold we used for the com- 
pactification should be S1/Z2 x Z'2 instead of a SijZ^- The behavior of the fields under the 
new Z'2 symmetry is, in fact, very similar to its behavior under the initial one. To see this 
let us define z' = z — Lq/2 and: 

z') = Le/2 + z!) = z) (39) 

Again, from Eq. (|l^), we can see that the Lagrangian is invariant under the Z2 symmetry: 

^(x", z') — > z') = r,t^(x", Le - z') (40) 

We will impose the same boundary condition as for Z2: 

V^(x°,;z') = ^(a;",^6 + -2') (41) 
Combining Eqs (B0[) and (HI]) we get: 



i)±{x'',-z') = ij^{x'^,z') 

ij±{x^,U-z') =i,^{x^,U + z') 



(42) 



which, in terms of ip and z, become: 



i}±{x°',z) = ij^{x°',Lfi - z) 
i)±{x'^,-z) = i;:f{x°',LQ + z) 



(43) 



Using this second parity we can find an explicit realization of the Q-symmetry as follows. 
First, we shall define the behavior of the fermions under this symmetry in the region I as, 

ij'iz) = Q^{z) = TM^) (44) 

now, Eq. (^) relates region I one and region II of the orbifold - as it should be since the 
physical space in a. S1/Z2 x Z'2 goes from to L/2 - so in order for Q to be a symmetry of 
the Lagrangian the fermions have to satisfy, 

Qi){U -z) = i)'{u -z)= r,^'(z) = r,r7^(^) = -v^i){u - z), (45) 
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where in the second and last equahties we have used Eq. (|43|) which can also be written as 
z) = T^tpi^x", Lq- z). 

Notice that this realization of the Q-symmetry is only possible in an even number of space- 
time dimensions since it is only in this case that there exists a matrix which anticommutes 
with all of the gamma matrices of the algebra and which does not belong to the algebra. 

With the above definitions, it is straightforward to see that the Yukawa coupling 

CY = f mx'^, z)^{x'', U-z)- U - z)) z) , (46) 

is invariant under all Z2, Z'2 and Q symmetries where Q^ix'^^z^ = Furthermore 
the action of the three parities forbids the presence of another non-local hermitian term, 



i {iI){z)iI){Lq — z) + iP{Lq — z)iIj{z))^{x'^ , z). In fact, Eq. 43 renders the above term to be 
identical to zero. 

In terms of x aiid rj. Eq. (|46|) becomes 



^Yi = f iixix", z)r]{x°', Lq- z) - Lq - z)r]{x°', z)) ^^^^ 

+ {f]{x"', Z)xix", Lq - z) - f]{x", Lq - z)xix°' , z))}^{x°' , z) . 

As before, the minimum energy solution for $ is 

( $ ) = h{z) . (48) 
From (0) and (^TD, the equation of motion for the surviving zero mode C,x,o{z) has the form: 

-d.^xM + t fh{z)i^fl{LQ - z) = . (49) 
In order to solve Eq. (|49|) , we shall use Eq. (^) that, in terms of x and rj leads to, 

i^,o{L,-z)=i^fl{z). (50) 

Because of the factor i in Eq. (|49|) S,xfi{z) will not be localized along z. 
The solution to ( ^Q]) with the ansatz (|50| ) is now given by 

ex,o(^) = ^e-(^), (51) 

where 

s{z) = f [ dz'h{z'). (52) 
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Making the SHO approximation as used in the five dimensional statement to be 

justified below, the properly normalized wave function for C,x,o{z) would be 

Cx^oiz) = ^e^^'^' . (53) 



From the above solution for the zero mode in the 6th dimension, Eqs. ( |5T| , p3[ ) , we 
notice a marked difference with the 5-dimensional case: the zero mode wave function is now 
oscillating inside the thick brane, along the sixth dimension, while in the five dimensional 
case, its counterpart has a localized form along the fifth dimension. 

Let us assume there is a kink solution for $, i.e. 

h{z) = V tanh(/i2;) , (54) 

where fi = (A/2)^/^f. With this solution ( p4D put into (p2D, the explicit expression for the 
non-vanishing zero mode is now 

V-^6 

Just as we have done with the five dimensional case, one could generalize the above 
discussion to include a "mass term" so that / h{z) —* f h{z) —m. As a result, one now has 

V-^6 

This more general expression (|56|) in fact determines the phase of the oscillation. 

In the construction of the mass matrices in four dimensions, we will need overlaps of 
wave functions in the extra dimensions, as we have discussed above in regards with the fifth 
dimension. How the mass matrices look like in six dimensions is the topic which will be 
discussed next. 

We end this section by presenting another type of Yukawa coupling which is used to 
actually localize fermions along the fifth dimension. The only difference with the previous 
section is that we now write it using the full six dimensions. With Tj defined in Appendix 
A, the appropriate coupling is 

SYuk2 = j rfVVrr^V- (57) 

Defining 75 = iVyV-j, one can see that Eq. (|57|) is invariant under ip{x^^y^z) — >■ 
±75^/'(a;'^, L5 — z) and $'(a;^, z) — $'(a;'^, L5 — z) which finally gives iIj±{x, —y, z) = 
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^Ibi'Ti^yVy ^) s^nd —y,z) = —^'{x,y,z). Also Eq. (|57D is invariant under the Q sym- 
metry provided that Q ^'{x^, y, z) = — y, z). Notice that Eq. ( pTj) can also be written 
as 

SYuk2 = j d'xf'ix'^'x - r]<^'v) ■ (58) 
Eq. ( p8|) will reduce to the usual coupling in five dimensions. One last commentis in 
order. Eq. ( ^Tf ) is also invariant under a simultaneous Z2-transformation: iIj{x",z) — *■ 
Tzi'{x°', Lq — z), $'(x, y, z) — > $'(x, Le — -z), as well as under the Q-symmetry. 

Before leaving this section, we would like to make a remark concerning Eq. (^31). Basically, 
it is a "mapping" of region I into region II and vice versa, namely iIj{x°', z) = T^iplx'^, Lq — z) 
or iP{x"',Lq — z) = V z%Ij{x" , z). Now, let us remember that Eq. (^31) is a consequence of 
our boundary conditions. When we substitute it into Eq. ( |^ ) so that one deals with the 
physical space which is now ranging from to it acquires a Lorentz non-invariant 

form ?/'(x", z)r^-?/'(a;", 2;). What this says is that our boundary conditions break the six- 
dimensional Lorentz invariance down to a five-dimensional Lorentz invariance. Our original 
Lagrangian (^Bj) is Lorentz invariant under the full six-dimensional Lorentz group and only 
when one goes to the physical space dictated by the boundary conditions, the six- dimensional 
Lorentz invariance is broken down to the five-dimensional one. 



B. (Almost) Pure Phase Mass Matrices 

We shall use the same notations as in Section ( |I B| ) . The action for the Yukawa interaction, 
in six dimensions, between the quarks and the SM Higgs field, is written as (the Down sector 
is treated in exactly the same manner) 

where Cg = ror2r2. We have, for the moment, omitted to write down other possible terms 
which are needed to determine the phases along the sixth dimension. This will be dealt with 
in the next section. We first begin with a "phenomenological" analysis. 

The previous analysis led us to write a generic (zero-mode) fermion field as 

^{x,y,z)=^{x)Uy)Uz), (60) 

Before making use of Eq. (^) to construct the mass matrix, let us describe a possible 
"geography" of the fermions along the extra dimensions. The discussion of Section ([[ B|) 
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pointed out the following features: The localization, along the fifth dimension y, of Qi at 
one place and f/f at another place produces a Democratic Mass Matrix as shown in Eq. (p!6[). 
That is the "geography" along the fifth dimension that we would like to keep. Basically, 
left-handed and right-handed fields are localized by two domain walls at different locations. 
Why this should be so is beyond the scope of this paper. However, one important point that 
should be kept in mind is the fact that, in our model, there are only two locations (left and 
right) along the fifth dimension, regardless of the family index, for each quark sector (Up or 
Down). As mentioned above, this gives rise to the universal effective Yukawa couplings gy^u 
and gY,d which determine the overall mass strength for each sector. Let us recall that gY,u 
and gY4 are proportional to the overlap between left and right for the Up and Down sectors 
respectively. Again, what splits gY,u from gy^d is beyond the scope of this paper. However, 
we will make some remarks concerning this issue at the end of the paper. 

The next question concerns the locations of various domain walls along the sixth dimen- 
sion. At the end of this section, we will present a simple example which shows how one can 
localize these domain walls. For the moment, we will simply parametrize these locations as 
shown in Eq. (0). We will assume that the domain walls which "fix" the phases for the 
three families are located at different positions along z. For the purpose of illustration, we 
will stay with this simple picture of family breaking in this manuscript. A more general case 
with phenomenological applications will be dealt with elsewhere. This will involve different 
profiles for different family kinks, etc.. 

We shall discuss below the implications of the cases when, for each family, Q and are 
"in phase" and when they are slightly "out of phase". But, first, let us use Eq. (|^) and 
Eq. (|59D to construct a general generic mass matrix for the Up sector. The mass matrix for 
the Down sector will be obtained in exactly the same manner. 

In the following, the quantity Lq which appears in various formulas is a generic symbol 
for the length of the physical space, which is Lq itself for the orbifold S1/Z2 or Lq/2 for the 
orbifold Si/{Z2 x Z'^). 

To begin, we will assume the following situation for the "geography" of family domain 
walls along the sixth dimension z. We will then discuss special cases of such a scenario. (As 
we have briefly mentioned above, this scenario is presented for the purpose of illustration 
and is not the most general case.) Let us define the following quantities which appear in 
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Eq. m-. 



(61) 



where i = 1,2,3 denotes the family index and where /ij = (A/2)^/^fj. Notice that, in prin- 
ciple, the quartic coupling A can depend on the family index i. This more general case, 
however, will be investigated elsewhere. From Eqs. one can write an effective 

Yukawa interaction in four dimensions and construct a mass matrix as we had done earlier. 
This construction is identical to the five-dimensional case, except that now the matrix ele- 
ments will contain an extra factor which is the overlaps of S,&{zys. As usual, the mass matrix 
will be similar to Eq. ([16D except that now, instead of the matrix elements being unity, one 
has 



where 



^33 



^ On ai2 ai3 ^ 



0-21 0-22 ^23 
0-31 0-32 033 



(62) 



(63a) 



-— / dz exp(i(mj+ — mjj)z) 
Le Jo 

{expiiinij^ —mj_)LQ) — l)/i{mj+ 



mj_)L6 , 



dz ^6,i+^6 ,i 
La 



(63b) 



1 f'^'' 

= — dz exp{i{ajln{cosh{^jz)) — ailn{cosh.{fiiz)) + {mi+ — mj )z) . 

Notice that Lg here is a generic symbol for the length of the physical space as we have 
mentioned above. 

The above equations (^ |63a| , |63b| ) refer to the case where domain walls, which "deter- 
mine" the phases of the fermions, are "located" at different places. We will specialize below 
to a few interesting possibilities. However, some important remarks can already be made. 
We ask the following question: Under what conditions will the mass matrix be hermitian or 
non-hermitian? 



1. Hermitian and Non-Hermitian mass matrices 
We now present two different scenarios. 
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(a) The parameters mi± which determine the "locations" of the domain walls possess 
interesting features. The first observation one can make is as follows. If the domain walls 
which "localize" the phases of Q and Lf^ (Left and Right), for each family, are located at 
the same place along z, i.e. 

mi+ = mi_ (64) 

one obtains the following results 

ttjj = 1 , (65a) 



ttji = a*j . (65b) 

The mass matrix Ai is hermitianl The hermiticity of the mass matrix is a consequence 
of the "collapse" of left and right (or Q and If^) , for each family, into the "same position" 
along the sixth dimension. Two remarks can be made concerning a hermitian matrix. First, 
its determinant is real. This means that arg(det A^) = 0. The possible connection of this 
statement with the strong CP problem (see e.g. a review by f^) will be explored further at 
the end of the paper. 

Let us first see if the hermitian matrix above is of a pure phase form. 

The discussion which follows will deal with issues which are also relevant to the non- 
hermitian case. 

Let us look at 

1 f^" 

ttij = — dz exp(z(aj ln(cosh(/ij2;)) — ln(cosh(/ij2;)) + (mj — 'mj)z)) . (66) 
-^6 Jo 

Under what conditions would a^j's look like pure phases, namely of the form e*^, or an almost 
pure phase of the form (1 — p)e*^ with p <^ 1? To answer this question, let us make a little 
detour to the meaning of wave function overlaps, thickness of domain walls and size of the 
extra dimensions. 

We have seen how one can localize fermions along the fifth dimension (y) by having do- 
main walls of sizes l/p -C L^. The effective strengths of various interactions are determined 
by the overlaps of the wave functions along y. For this reason, it is preferable to have the 
thickness of the domain walls small enough, i.e. 1/p ^ L5, so one can "fit" several fermions 
along y in such a way as to obtain desirable effects such as "slow" (or no) proton decay, 
possible mass hierarchies between different fermion sectors (quarks, leptons), etc... As we 
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move on to the sixth dimension, it is not obvious that such a picture is still necessary. In 
fact, at least as far as the pure phase mass matrix is concerned, the thickness of these domain 
walls can be as large as the size of the compactified dimension itself, as we shall see below. 

Let us, for the time being, assume that all domain wall thicknesses (along z) are of the 
size of the compact dimension, i.e. ~ 0{Lq). In this situation, one can use the SHO 
approximation and carry out the integration of Eq. (pB]), namely 

Qij = — / dz exp(— z(A/i?-2;^ — Arriijz)) , (67) 
Lq Jo 

where 

^f^lj = (l/2)(a^M• - dj/if) , (68a) 



Arriij = rrii — rrij . (68b) 
The integration can be explicitly carried out. One obtains 

I ^(2A/ifj.L6 - Aniij) \ ^^^^i i Ami,- 



= — == exp i . (69) 



In a phenomenological application of Eq. (^9]), one can use it without making any approx- 
imation. However, in order to see if it has a more familiar pure phase form or not, we will 
make an expansion of (|69D . 
Let us define 



A/i?. Le = X,, , (70) 



Arriij Lq = Vij . (71) 
For Xij, Uij < 1, one can expand ( |5^ giving 

where we have neglected terms of 0{xfj,y'^) or less in the modulus and terms of 0{xfj,y'^) 
in the phase. Notice that for aji, one has x|j = — x'^j and yji = — yij, and hence a^j = a*j as 
they should. In this form one can see that the hermitian mass matrix is almost of the pure 
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phase form. This would have been the case if one could neglect terms containing Xij and i/ij 
inside the coefficient multiplying the exponential. However, we will not neglect those terms, 
leaving the possibility of a small deviation [^] from a pure phase mass matrix. 

Notice that when the domain walls are all located at the same point, i.e. Arriij = 0, and 
when they have the same thickness (or /ij = fij), one recovers the DMM form, namely 
aij = 1, as one can see from Eq. |^J7^ ). In addition, we notice that one can also obtain 
the almost-pure phase hermitian mass matrix when either Arriij 7^ or Afifj 7^ 0, but not 
necessarily both, as can easily be seen. 

(b) As we have seen above, within the framework of Eqs. ( |63a| , |63bD , the mass matrix can 



be purely hermitian provided the condition mj+ = mj_ is fulfilled. What would happen if 
mj+ 7^ To study this question, let us refer back to Eqs. ( |634|63b| , [f^ ) and let 



mi+ - mi_ = ei . (73) 



Also for convenience, let us define 



5i = eiLe. (74) 

With the above definitions, the diagonal matrix elements which are no longer unity, can be 
written as 

an = exp{tSi/2) . (75) 

{di/2) 



The off-diagonal elements are similar to Eq. (|7^), except that now one has the following 
replacement y-ij Ui+j- = (^i+ — nijJ)L^. It is convenient to remove the phases from the 
diagonal elements by absorbing the phases into ^6,i+5 namely ^6,i+ = exp(z(5j/2) ^6,i+- From 



the definitions of an and a^j, one now has 



sm(^./2) 



2 1 

^ 45 24^*-'^' 



1 

12 



^ijUi-J- 



1 A 1 

12 24 



5f + ^45.}exp0(^-^)}, 

(77) 



where 



Notice that Vi-j- 



Vi-j- = {mi_ - mj^)LQ . (78) 
-yj-,i-- In Eq. (|77D , we have made use of the above phase redefinition 



and of (|73| , [74D . The mass matrix described by the above elements is not hermitian for 
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the following reason. The modulus of aji will have a term —j^yj^^i^Sj — + j2^^i^j = 
— — ^a;fj-5j. It can easily be seen that \aij\ ^ |ajj| unless 5j = — 5j which cannot 
be satisfied for all j. Despite the fact that the phase of a^i is the negative of that of a^^ 
the difference in modulii implies that, in general, a^j 7^ a*j-, and hence the non-hermiticity 
of the matrix. It can be approximately hermitian if one can neglect the terms containing 6i 
in (0). 

Notice that, even for the special case where all "left-handed" family domain walls are 
"located" at one point along z, i.e. mj_ = mj_ = m_, so that Vi-j- = 0, and all "right- 
handed" family domain walls at another place, i.e. 6i = 6, the non-hermiticity still appears 
in the difference in mudulii between aji and aij because of the presence of 6. 

From the above discussion, one can see that one recovers the hermitian matrix in the 
limit 6i 0. 

In summary, we have shown that, in general, the deviation from hermiticity in our frame- 
work comes from the splitting between "left" and "right" , namely mj+ 7^ mj_ . 

The above analysis can be carried over to the Down sector in exactly the same manner. 
There are however two interesting remarks that can be made. First, although the mass 
matrix for the Down sector is now characterized by a universal strength gy^d which is in 
general different from gy^u, the matrix itself can be identical to the one for the Up sector if 
we consider scenario (a). The reason is that scenario (a) is one in which the domain walls 
for Q and V^, for each family, are "located" at the same place along the sixth dimension, 
which is exactly the same as for the Up sector. Therefore the matrix elements (without the 
universal strength) are the same. In consequence, the diagonalization matrices are the same, 
i.e. Vu = Vd. Hence, Vckm = Vu^d = 1, a mere unit matrix. In other words, the mass 
matrices for the Up and Down sectors cannot be both hermitian. To obtain a non-trivial 
CKM matrix, at least one of the two matrices has to be non-hermitian in this particular 
scenario. 

The above (almost) pure phase mass matrix as obtained from six dimensions is what we 
have set out to derive. From it, we have learned a few things. 

(a) In general, the almost pure phase form of the mass matrix can be easily seen if the 
thickness of various domain walls along the sixth dimension is of the order of the compact- 
ified sixth dimension. (There is no reason why, in principle, the thickness of the domain 
walls should be much smaller than the compactified dimension, in contrast with the five- 
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dimensional case.) 

(b) When the domain walls "fixing" the phases for Q and Lf^, for each family, are located 
at the same place, (m„ j+ = mu,i-), the mass matrix is purely hermitian. As we have seen 
above, another possibility is when the domain walls "fixing" the phases for Q are at one 
location and those which are responsible for "fixing" the phases of If^ are at another location, 
in which case the mass matrix is also hermitian. If one considers these cases to be a "tree- 
level" situation -a statement to be further clarified below, the fact that arg(det Al) = 
makes this scenario an interesting "candidate" for a solution to the strong CP problem. 

(c) The mass matrix becomes non-hermitian when mu,i+ 7^ iriu^i-- We will briefiy discuss 
below the possibility that mu,i+ 7^ is due to "radiative corrections" of the case m„^i+ = 

The mass matrix for the Down sector is obtained in a similar way. The main difference 
between the two sectors is the "universal" strength which appears in fron of the matrix: 
9y,u-^ for the Up sector and gY,d-^ for the Down sector. The other difference in the case 
of a non-hermitian matrix (Scenario (b)) is the sphtting between "left" and "right" for each 
family, which does not have to be the same for the two sectors. 

Notice that, in order to be more general, wc aloow the possibility of two different mass 
scales: and v^- If there were only one SM Higgs field then Vu — Vd — v. In this case, the 
disparity between the mass scales of the Up and Down sectors would come from from the 
difference between gy^u and gY,d, which, in turns, could come from the differences between 
wave function overlaps, along the fifth dimension, of the two sectors (modulo differences in 
the fundamental Yukawa couplings). To keep our discussions as general as possible, we also 
allow for the possibility that two SM Higgs fields exist. 

It is beyond the scope of this paper to discuss in detail the phenomenology of our model. 
It will be carried out elsewhere. 

2. Some Remarks on localization of family domain walls along the sixth dimension 

In this section, we will briefiy discuss one way to localize the various domain walls re- 
sponsible for "fixing" the phases of fcrmions along the sixth dimension. There are probably 
several mechanisms to achieve this. We will present one of such mechanisms, from the point 
of view of effective field theory. 
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For simplicity, we shall assume in this section that = a = f/^/\/2. This simple 
assumption basically refers to couplings between fermions and background scalar fields which 
are invariant under the family symmetry. 

First, let us list the parameters that we need to construct an almost pure phase mass 
matrix. From Section ( [II B 1|) , we learned that we need: /ij with i = 1, 2, 3 which control the 



thicknesses of the domain walls and mi± which control the locations of the domain walls. 
We also learned that, one can obtain a hermitian mass matrix when mj_|_ = mj_ and a 
non-hermitian matrix when mj_|_ ^ rrii-. It turns out to be a highly non-trivial task to find 
a mechanism which can "explain" the origin of these parameters. In some sense, it might 
even be overly ambitious to make such a claim. We will, however, make an attempt to, at 
least, hint at one possible scenario. 

In Section (|IIB 1|) , we were basically doing the "geography" of family domain walls along 



the sixth dimension. To construct a scenario for the "geographical points" (the various m's), 
let us recall that the family symmetry of our model is 5*^ ® S^" . The background scalar 
fields which couple to Q or If^ will appear in terms such as Q^Q, U'^^U'^. We will therefore 
need two of such background fields in order to write down invariant Yukawa couplings: $q 
and $c/c. These background fields, and $[7c, will be represented by 3 x 3 matrices. Some 
of the details concerning the potential for these scalars are given in Appendix |B[ Here, we 
will just quote the results. The discussion below refers to the Up sector. As we have seen 
earlier, the Down sector can be treated in exactly the same manner. 

We will concentrate on scenario (a) of Section ( [II B 1| ) for the purpose of illustration. We 
will assume the following Yukawa interactions: 

J^Y = fQ^qQ + fU'^u^U' + h.c. , (79) 

where, for simplicity, we have put the two Yukawa couplings to be equal. (A more general 
case can be accommodated straightforwardly.) The minimization of the potential gives, at 
tree levels 



( hi{z) ^ 



(80) 



h2iz) 
\ h{z) J 

One could assume that, at some deeper level and because of the family symmetry, the two 
background fields behave in exactly the same manner, i.e. having similar parameters, and , 
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in consequence, one has 



( hi{z) ^ 



$c/0= h2{z) . (81) 
^ hiz) ^ 

These VEV's will be shifted by radiative corrections. It is beyond the scope of this paper 
to examine this problem and we will simply parametrize these shifts by 

h,{z) hi{z) + 5h, , (82) 

where the shifts are assumed to be independent of z and are also assumed to be much smaller 
than Vi (or /ij). 

Combining Eq. (^) with Eq. (^), one can make the following identification 

rrii- = mi+ = fShi . (83) 

This is the case when one would obtain a hermitian mass matrix of scenario (a) of Section 
( [11 B 1| ) ! It goes without saying that there are two assumptions which have been made. First, 



we have assumed the equality of the Yukawa couplings in Eq. (^). Second, we have assumed 
that the behavior of the two background scalar fields are identical. These assumptions might 
come from some deeper symmetry between Q and If^ (or Z)^). This is very similar to the 
notion of left-right symmetry that one encounters in four-dimensional model building. In 
consequence, the hermiticity of the mass matrix that we obtained by "phenomenologically" 
putting mj_ = mj+ might be justified by some form of left-right symmetry. 

In addition (|82D , one should also take into account vertex corrections which will be differ- 
ent for Q and (they have different gauge interactions for example). Let us parametrize 
those shifts by 

fQ = f + 5fQ, (84a) 

fu = f + Sfu^ , (84b) 
where the notations are self-explanatory. We will assume that Sfq^w^ -C /. Naturally, 
/q 7^ fu- 

From the above equations, one can make the following identifications: 

mi_ = fQ6hi , (85a) 
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mi+ = fuShi (85b) 

Since one expects Sfq ^ 6fu<^ and, in consequence, /g 7^ fu, one would expect, in general, 
mj+ 7^ mj_ which is a condition for the appearance of a non-hermitian matrix. However, an 
approximate hermitian matrix could arise if the radiative corrections and, in particular, the 
difference in the radiative corrections are small. One can see that, as we turn off whatever 
interactions (gauge, etc..) which contribute to the vertex corrections Sfg^u", one recovers 
the hermitian case, namely mj_ = mj_|_. 

Pursuing the same idea, one can also assume that Z^^'s have a similar coupling of the 
form fD'^f^ccD^. Assuming that {^d'^) has a similar form to Eqs. ( pO| , |8TD , one can now 
see that, in the absence of vertex corrections, one obtains mj_ = mj.„^+ = mi^d,+, which is 
just scenario (a) discussed above. Since -D^ and have different quantum numbers, one 
expects that their vertex corrections will be different from each other. In consequence, one 
will obtain mass matrices of the form (32) with coefficients of the form ( |63a , p3b ). 

In the scenario just outlined above, one can make interesting connections with the strong 
CP problem. In the absence of vertex corrections, the mass matrix is hermitian and hence 
arg(det7W) = 0, a possible solution to the strong CP |§] problem? (One could assume 
CP to be a symmetry of the Lagrangian so that Oqcd = 0.) As mentioned above, this 
hermiticity might come from some left-right symmetry {Q ^ W^, D'^) which gives mj_ = mj+ 
at "tree level". It could be quite provocative to see if there are connections, if any, with 
previous solutions to the strong CP problem which made use of the quintessential Left-Right 
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Turning on the vertex corrections, the pure phase mass matrix becomes non-hermitian 
and, as a consequence, one would obtain a non-zero contribution to the strong CP parameter 
6. If this were truly a plausible scenario for the strong CP problem, the resultant 6 should 
obey the upper bound of ~ 10^^. However, it is beyond the scope of this paper to analyze 
its magnitude. We will come back to this issue in a subsequent paper. Our future studies 
will focus on the following two questions. Will the "radiative corrections" be small enough 
so as to account for both the phenomenological constraints on the mass matrices and the 
magnitude of 67 If those phenomenological constraints on the mass matrices require a 
"large" radiative correction, is there a "natural" mechanism to make 6 small enough? 
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III. CONCLUSION 



In this paper, we have studied the problem of fermion mass hierarchy from the point 
of view of large extra dimensions. To this end, we have added two extra compact spatial 
dimensions. In particular, we have shown how one can construct a particular kind of mass 
matrices which is very successful in fitting the pattern of quark masses and mixing angles: 
The pure phase mass matrix. This matrix is characterized by a universal Yukawa strength 
appearing in front of a matrix whose elements are of the form exp{i6ij). In our construction, 
the universal Yukawa strength arises from the overlap of the wave functions of the left- 
handed quarks (denoted by Q) and the right-handed quarks (denoted by and D'^) along 
the fifth spatial dimension (y). Along y, all left-handed families are localized at one place 
and all right-handed families at another place, with the localization carried out by domain 
walls whose thicknesses are assumed to be much smaller than the radius of compactification 
of y. We then proceed to show that the almost pure phase mass matrix arise from the overlap 
of wave functions between different families and also between left-handed and right-handed 
quarks, along the sixth dimension z. Along z, the "phase determination" is carried out by 
domain walls whose thicknesses are assumed to be of the size of the radius of compactification 
of z. 

The almost-pure phase mass matrices obtained in six dimensions have some interesting 
properties, according to the "locations" of the family domain walls, which fix the phases, 
along the sixth dimension. In one case (scenario (a)) which is dubbed "tree level" in this 
paper, the domain walls for Q and If^ or D'^ are "located" at the same place along the 
sixth dimension z, for each family. The mass matrices thus obtained are purely hermitian. 
In addition, apart from a different universal Yukawa strength, the matrices of the Up and 
Down sectors are identical, giving rise to a situation in which the CKM matrix is simply a 
unit matrix. We then considered a scenario in which either the domain wall for If^ or D'^, or 
both, is split from that for Q. As we have shown in Section ( [II B 1|) , this would imply that 
the mass matrix of at least one of the two sectors is non-hermit ian, and the two matrices 
will be different from each other, implying a non-trivial CKM matrix. 

One should also keep in mind the possibility that the mass matrices of both sectors are 
hermitian but not identical. In this type of scenario, one would get a non-trivial CKM 
matrix as well as a correct spectrum. This possibility is mentioned at the end of Section 



30 



( [II B 1| ) . In order to be able to build a model for the "locations" of various quarks along the 
extra spatial dimensions, a full phenomenological analysis of various possibilities should be 
carried out to serve as a guidance. 

These two cases of hermitian and non-hermitian mass matrices might have important 
connections to the strong CP problem as we have briefly discussed above. This interesting 
issue will be further investigated in a future paper. 

Finally, a number of interesting issues such as the Kaluza-KIein modes, the extension to 
the lepton sector, and others will be dealt with in future publications. 
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APPENDIX A 



In this appendix we are going to present a brief review of the spinorial representations of 
the orthogonal group, 0{D), in higher dimensions {D > 4). We are going to follow closely 



the treatment done by Weinberg in his book [14|, with a slightly different notation. 



The starting point is a set of matrices, which spawn the Clifford algebra, with the anti- 
commutation relations {7^,7;/} = 25^,^. In addition, our attention will be fixed on spaces of 
even dimensionality {D = 2n) and, subsequently, we'll extend it to odd dimension spaces. 

Using the anticommutation relations of the 7 matrices, we can define n fermionic harmonic 
oscillators a.s af = ^ {—l2i + ^72i+i) with i = 0, . . . , n~l that are independent and, therefore, 
the set of basis vectors of the representation space has 2" elements which can be written as: 

|siS2 ...Sn) = a^^'a^'"' ■ ■ ■ a+''" |0) (Al) 
being |0) a vacuum annihilated by all destruction operators Oj. In this basis the matrices 
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take the form: 

/-I o\ [-1 o\ [o l\ 

ai= ® ® #2x2 • • • ® i^2x2 (A2) 

\o l) \o l) [(,(,) 

the — I's are due to the fact that and anticommute. Finally the 7 matrices can be 
easily obtained and they read as: 

72i = -C^S <8) • • • <8) CTg (g) CTi (g) #2x2 <8) • • ' <8) -^2x2 . 
72i+l = (T3 ® ■ ■ ■ ® 0-3 ® (72 ® -^2x2 ® • • • <8) l'2x2 

where CTj's are the Pauli matrices. Note that this representation does not give the usual 
representation in four dimensions but we can relate both using the following unitary trans- 
formation: 

-^-^((73 + (72)®---® ((73 + (72) (A4) 

SO the "usual" representation is 

72i = (^2 • • • (72 <8) CTi (8) l'2x2 ' ' ' <8) l'2x2 , , , 

(A5) 

72i+l = (72 (8) • • • (8) (72 (8) (73 (8) l'2x2 <8) • ' • <8) l'2x2 

Since, 

72i72i+l = i -^2x2 <8) • • • ® -^2x2 ® Cr2 ® -^2x2 ® " " " ® -^2x2 (A6) 

the product of the 2n 7 matrices is: 

2n-l 

n 7i = ^"f^2®---®(T2 = ?772n- (A7) 

Where 1] is a phase such that 72„72n is the identity. Note that we are labeling the gamma 
matrix equivalent to 75 with 2n instead 2n + 1. This difference comes from our choice for 
the labeling starting from instead of 1. This new matrix anticommutes with all 7's and 
therefore it implies that all spinorial representations of 0{2n) are reducible. 

Let us find now the spinorial representations of the orthogonal groups with odd dimen- 
sionality, D = 2n+ 1; this is much more simpler once we have the representation for 0{2n) 
since we just have to take this representation and add the 72^ matrix. In this case the repre- 
sentation is irreducible because we can not find any independent matrix that anticommutes 
with all the gamma matrices. 

The transition of the 0{D) representations to — 1, 1) representations is done through 
a wick rotation. 

To finalize we will explicitly write the gamma matrices for 0(5, 1). 
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Six dimensions (with metric ( — h + + ++)): 

To = i (72 <8) (7i (8) 1^2x2 = 



r2 = (J2 (8) (J2 (8) (72 



= (72 (g) (73 (g) #2x2 



04x4 -70 
70 04x4 

04x4 —72 

72 04x4 

'^04x4 -ilb 
yijb 04x4 

= (T3 (g) #2x2 ® #2x2 




'4x4 


-71 


71 


04x4 


4x4 


-73 


73 


04x4 


■( 


04x4 




#4x4 



-^4x4 04x4 
04x4 —#4x4 



(A8) 



Notice that, in the above equations, 7^* {/i — 0,1, 2, 3) and 75 are simply defined here as 
7o = (7i (g) #2x2, 7i = «C2 <g) Cj, 75 = (73 (g) #2x2- Thcse definitions just happen to coincide 

with the 4-dimensional ones with a metric (H ). This is simply a compact way of 

writing the 6-dimensional F's. There is no change in metric. To see how the 6-dimensional 
metric ( — h + + ++) reduces to a 4-dimensional metric ( — h +- 1-) when the two extra 
spatial dimensions are compactified, one rewrites 7'* in terms of the gamma matrices which 
correspond to the metric ( — namely 7^ = ^7^. In this way, the kinetic terms will be 
preceded with a plus sign when they are reexpressed in terms of 7'^. Of course, 75 remains 
unchanged. 



APPENDIX B 



In this article, we had been discussing models in which fermions are localized at one 
place or another along the extra dimensions and inside fat branes with the same or different 
widths and how these settings could affect the phenomenology of the 4D models. However, 
we did not provide any model that explains these different settings; this will be addressed 
in this appendix. 

As an example we are going to study the possibility that the background scalar field is a 
composite of fields that transforms under a three dimensional representation of the family 
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group; therefore this background scalar field $ takes the form: 



( 01 (x^^ 



= (01 (x) 02(X) 03(X))+ ® 



(Bl) 



02 W 

\03(a;) y 

where the 0j's are the "fundamental fields" from which the background scalar field is com- 
posed. 

The first of the models we are going to propose consists on a 0"^ potential, without cubic 
terms, for 2 composite fields of the form (|BI|), 

y($i,$2) = ^Tr[$i$+] + ^Tr[$2$^] + ^(Tr[$i$+] + h.c.) 

+ ^ Tr[($i$+)2] + ^ Tr[(<|.2$+)^] + ^ Tr[($i$+)($2$+)] 



+ ^ (Tr[$i$+$i$+] + h.c.) 

2 2 

= ^l^il + + m3|<l>i||<l>2| cos^a 



(B2) 



Ai 



A, 



A, 



+ -i |<l>i|^ + |<l>2|^ + ^ |$l|1$2r COS^ « 

+ ^|$in$2|'cos^a 



where n\\ 2 are negative coefficients. We also made use of the following definitions. 



l$l 



Tr[$$- 



cos a = Re 



Tr[<l>i$5 



(B3) 
(B4) 



If 2m3 > — (A3 + A4 cos^ a)|$i| |$2| holds then the minimum of the potential occurs when 
both fields take expectation values along orthogonal directions, that is, when Tr[$i$^] = 0. 
In consequence, if we suppose that the expectation values for |$i| and |$2| are u and v 
respectively, the potential we have to minimize is: 



, , TTii o Tfln r) Ai /I At a 



(B5) 



and its minimum is located at: 



u 



7712 



(B6) 
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therefore we can suppose that, 



u 









V 










(B7) 



(B8) 



This means that this model will localize the components of the family multiplet at different 
positions along the sixth dimension, namely 0, u and v. 

We can extend this model to localize all the components of the family multiplet inside 
the orbifold, with the background scalar field in (^6[) having all three eigenvalues different 
from zero. This can be done if we add another background field that can be a singlet, or 
a composite like the ones used. In the former case the three components will be shifted by 
the same amount, s, ending in positions s + u, s + v and s; in the later case the fermions 
will be located at u, v and z (the expectation value for the third field). 
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